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Controllability and observability of grid graphs 
via reduction and symmetries 

Giuseppe Notarstefano Gianfranco Parlangeli 



Abstract 

In this paper we investigate the controllability and observability properties of a family of lin- 
ear dynamical systems, whose structure is induced by the Laplacian of a grid graph. This analysis 
is motivated by several applications in network control and estimation, quantum computation and 
discretization of partial differential equations. Specifically, we characterize the structure of the grid 
eigenvectors by means of suitable decompositions of the graph. For each eigenvalue, based on its 
multiplicity and on suitable symmetries of the corresponding eigenvectors, we provide necessary and 
sufficient conditions to characterize all and only the nodes from which the induced dynamical system 
is controllable (observable). We discuss the proposed criteria and show, through suitable examples, how 
^- such criteria reduce the complexity of the controllability (respectively observability) analysis of the grid. 
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^ I. Introduction 
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In several modern engineering application areas, there are important physical phenomena 
whose dynamic model is induced or strictly related to the structure of a graph that models the 
interaction among components of the main system. 

For example, in multi-agent system control (e.g., distributed robotics, sensor networks or 
smart power grids) a communication or interaction graph induces the structure of the feedback 
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control input. In Markov chains the evolution of the probabilities of a finite number of states 
can be modeled as a dynamical system structured according to the graph of one-step transition 
probabilities. Similarly, in quantum computation the evolution of interacting particles obeying to 
quantum laws may be described by differential equations structured according to an interaction 
graph. Another area in which a dynamic model is related to a graph structure is the one of 
discretized partial differential equations. In this case the graph determined by the discretization 
rule induces the structure of the approximating ordinary differential equations. 

In this paper we will concentrate on linear time invariant dynamical systems whose state 
matrix is induced by the Laplacian of a fixed undirected graph. In particular, we consider d- 
dimensional grid graphs, also know as lattices. This graph topology appears in many important 
application scenarios as we show in the next sections. 

We investigate the eigenstructure of the grid graph Laplacian in terms of structural properties 
of the graph, and study how this structure affects the controllability and observability properties 
of the induced system. Our main goal is to relate the controllability and observability properties 
to graph theoretic rules involving simple arithmetic operations on the graph labeling. 

Controllability of complex network systems has received a widespread attention in the last 
years in several areas [3]. We organize the relevant literature in three parts according to the three 
main motivating scenarios for our problem set-up. First, a system with Laplacian state dynamics 
arises in network systems running an average consensus algorithm. A survey on these algorithms 
and their performance may be found in 01 and references therein. The controllability problem for 
a leader-follower network was introduced in (51 for a single control node. Intensive simulations 
were provided showing that it is "unlikely" for a Laplacian based consensus network to be 
completely controllable. In [6j, see also flTj, "necessary and sufficient" conditions are provided 
to characterize the controllability and observability of path and cycle graphs in terms of simple 
rules from number theory. In [[8] and ||5]|, see also |[T0l , necessary conditions for controllability, 
based on suitable properties of the graph, are provided. Other contributions on the controllability 
of network systems can be found in ifTTTl . Ifl2l . lfT3~l . Observability has been studied for the first 
time in |fl4|. where necessary conditions for observability, as in the dual controllability setting 
investigated in flU and 0, are provided. A parallel research line investigates slightly different 
properties called structural controllability, [fT5l . QUI, and structural observability, [fTTll . Here, 
the objective is to choose the nonzero entries of the consensus matrix (i.e. the state matrix of 
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the resulting network system) in order to obtain observability from a given set of nodes. It is 
worth noting that controllability and observability of a network system are necessary structural 
properties in many interesting network problem as estimation, intrusion detection and formation 
control, e.g., 0, ttH, EH, HH- 

Second, continuous time quantum walks can be modeled as linear time invariant systems 
whose state matrix is the imaginary skew-Hermitian matrix iH , where H is called Hamiltonian 
and can be either the Adjacency matrix or the Laplacian of the transition graph fl2T|. 

The state transfer problem (which is strictly related to the controllability problem) for quantum 
systems is investigated in ||22|. The paper explores the eigenstructure of the Hamiltonian (which 
is taken as the Adjacency matrix of the underline graph) to characterize the state transfer. A key 
reference establishing a connection between our controllability analysis and the controllability of 
quantum walks is ll23l . Here the controllability of continuous time quantum walks is investigated 
and related to the controllability of a linear time invariant system with the structure considered 
in our paper. The controllability of continuous time quantum walks on graphs is also studied 
in ||24l| and ll25l . More specifically, <i-dimensional grid or lattice graphs play an important role 
in quantum computation. The controllability problem on this specific graph structure has been 
investigated in [|26l and ll22l . 

Third and final a system with the structure studied in the paper appears when discretizing 
partial differential equations (PDEs) containing the Laplace operator Il27ll . Such systems include 
several diffusion and wave propagation equations appearing in fluid-dynamics, mechanics, acous- 
tics and electromagnetism. In ll28l discretization of PDEs was indicated as a motivating example 
for the analysis of the Laplacian eigenstructure. The controllability of a discretized version of the 
heat equation on a one dimensional grid domain is investigated in [|29l and extended to the case 
of constrained input in [|30l . Finally, in OTl and ll32ll trajectory planning of multi- agent systems 
is performed by studying a partial differential equation describing a continuum of agents. That 
is, the multi-agent dynamics is obtained as a discretized version of a partial differential equation. 
Controllability is guaranteed by the particular choice of the control nodes. A more general choice 
of the control nodes leads to our controllability problem. 

The contribution of the paper is threefold. First, we identify a mathematical framework, namely 
the controllability and observability of linear time invariant systems induced by the Laplacian 
of a grid graph, that has numerous applications in several engineering areas. In particular, we 

February 24, 2013 DRAFT 



highlight how this framework appears in distributed control, quantum computation and discretized 
partial differential equations. 

Second, we characterize the structure of the Laplacian eigenvectors of a grid. Namely, we 
show that, on the basis of a prime number factorization of the grid dimensions, the eigenvector 
components present symmetries related to suitable partitions of the main grid into sub-grids that 
we call bricks. Given a partition of the grid graph into bricks, we show that the eigenvalues 
of the elementary brick are also eigenvalues of the main grid. Also, the grid eigenvectors 
associated to the common eigenvalues are obtained by composing (with suitable flip operations) 
the corresponding eigenvectors of the basic brick. Furthermore, we show that in each brick (and 
thus also in the main grid) the eigenvector components may show symmetries with respect to 
one or more of the grid axes. 

Third and final, we provide necessary and sufficient conditions to completely characterize the 
controllability and observability of grid graphs. We start showing that loss of controllability and 
observability can be studied by identifying all the zero components of an eigenvector. Based on 
the evaluation of suitable sets of polynomials, together with the eigenvector symmetries, we are 
able to determine all and only the eigenvector components that can be set to zero simultaneously. 
Thus, on the basis of the node labels, the eigenvector symmetries and the polynomial evaluations, 
we provide easily implementable routines to: (i) identify all and only the controllable (observable) 
nodes of the graph, (ii) say if the graph is controllable (observable) from a given set of nodes 
and (iii) construct a set of control (observation) nodes from which the graph is controllable 
(observable). 

The paper is organized as follows. In Section [XT] we introduce preliminary definitions and 
properties of undirected graphs, set up the controllability and observability problems and describe 



the motivating scenarios for our framework. In Section III we characterize the controllability and 
observability for grid graphs with simple eigenvalues. In Section [TV] we analyze the symmetries 
in the structure of the grid graph eigenvectors. On this basis, in Section [V] we provide necessary 
and sufficient conditions for the controllability (observability) of general grid graphs. Finally, in 
Appendix we recall results from [6] on the controllability (observability) of path graphs. 

Notation: Let N denote the natural numbers, for i e N we let e; be the i-th element of the 
canonical basis, e.g. e\ = [1 ... 0] T . For a vector v e R. d we denote (v)i the ith component 
of v so that v = [(f) i . . . (v)d] T - We denote II e R. dxd the permutation matrix reversing all the 
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components of v so that Uv = [(v) d . . . (v )i] T (the j-th column of n is [II] j = e n _ J+ i). Adopting 
the usual terminology of number theory, we say that k is a factor of m if there is an integer q 
such that m = kq. Given two integers b and c, if an integer m is a factor of 6 — c, we write 
6 m = m c. We denote GCD(a, b) the greatest common divisor of two positive integers a and b. 

II. Problem set-up and motivations 

In this section we present some preliminary terminology on graph theory, introduce the network 
model, set up the controllability and observability problems and provide some standard results 
for linear systems that will be useful to prove the main results of the paper. 

A. Preliminaries on graph theory 

Let G = (I, E) be a static undirected graph with set of nodes / = {1, . . . , n} and set of edges 
E C I x I. We denote Mi the set of neighbors of agent i, that is, Mi — {j E I \ (i,j) E E}, and 
dj = YljeAf- 1 me degree of node i. The degree matrix D of the graph G is the diagonal matrix 
defined as [D] u = dj. The adjacency matrix A E ]R nxn associated to the graph G is defined as 

14,-j 1 if< " J)e£ 

I otherwise. 

The Laplacian L of G is defined as L = D— A. The Laplacian is a symmetric positive semidefinite 
matrix with k eigenvalues in 0, where k is the number of connected components of G. If the 
graph is connected the eigenvector associated to the eigenvalue is the vector 1 = [1 ... 1] T . 

Next, we introduce the notion of cartesian product of graphs. Let G = (I, E) and G' = (/', E') 
be two undirected graphs. The cartesian product G □ G is a graph with vertex set I x I' (i.e. the 
cartesian product of the two vertex sets) and edge set defined as follows. Nodes [i,i'] E I x I' 
and [k,k'] E I x I' are adjacent in GOG if either i = k and (i',k') E E' or %' = k! and 
(i, k) E E. The cartesian product is commutative and associative. Thus, a d E N dimensional 
product graph, Yl e=1 Ge, is constructed by combining the above definition with the associative 
property. 

We introduce the special graphs that will be of interest in the rest of the paper. A path graph 
is a graph in which there are only nodes of degree two except for two nodes of degree one. 
The nodes of degree one are called external nodes, while the other are called internal nodes. 
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From now on, without loss of generality, we will label the external nodes with 1 and n, and the 
internal nodes so that the edge set is E = {(i, i + 1) | i 6 {1, . . . ,n — 1}}. 

A c/-dimensional grid graph (or lattice graph) is the cartesian product of d paths (of possibly 
different length). In a grid graphs the nodes have degree from d up to 2c?. We call the nodes 
with degree d corner nodes. Corner nodes are obtained from the product of external nodes in 
the paths. 

Given a rf-dimensional grid graph G = P\\3 . . . nP d , we denote i = [(i)i, . . . , (i)d] a node of 
G, where the component (i) K identifies the position of the node on the Kth path. Also, given 
a Laplacian eigenvector of the G, w E R m ~ nd , we say "the component [(i)x, . . . , (i)d] of w" 
meaning "the component (i) 1 ■ (rii ■ n 2 ■ . . . ■ n d ) + (i) 2 ■ (n 2 • ■ ■ ■ • n d ) . . . + (i) d of w". 

B. Controllability and observability of graph induced systems: problem set-up and analysis tools 

Next, we introduce the class of systems that we investigate in the paper. Informally, we 
consider linear time invariant systems whose state matrix is the Laplacian of a grid graph, the 
input matrix is obtained by directly controlling a subset of the node dynamics and the output 
matrix by observing a subset of the node states. Formally, let G = (I = {1, . . . , n}, E) be a grid 
graph, I c = {ii, . . . ,i m } C I and I = {ji, . . . , j p } C /, a first order dynamical system induced 
by G, I c and I is the system 

x(t) = aLx{t) + Bu(t), 

(1) 
y(t) = Cx(t), 



' -i 



andC 



Cjj 



iT 



fl h> 



where a is a scalar, L is the Laplacian of G, B 

It is a well known result in linear systems theory that the observability properties of the state- 
output pair (L, C) correspond to the controllability properties of the state-input pair (L T , C T ) = 
(L,C T ). Thus, the controllability and observability analysis for the class of systems in ([T]) can 
be performed by using the same tools. 

We start with some notation. The set of states that are controllable is the controllable subspace 
and will be denoted X c . Respectively, the set of initial states that produce an identically zero 
output is the unobservable subspace and will be denoted X no . 

An important result on the controllability (observability) of time-invariant linear systems is 
the Popov-Belevitch-Hautus (PBH) lemma, e.g. [|33l . Combining the PBH lemma with the fact 
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that the state matrix is symmetric (therefore diagonalizable) the following lemma follows. 

Lemma 2.1 (PBH lemma for symmetric matrices): Let A E R nxn , B G R nxm and C G R pxn , 
n,m,p G N, be the state, input and output matrices of a linear time-invariant system, where A 
is symmetric. Then, the unobservable subspace X no associated to the pair (A, C) (respectively 
the orthogonal complement to the controllable subspace X c associated to the pair (A, B)) is 
spanned by vectors vi satisfying for some AgR 

B t vi =0 Cvi =0 

respectively (2) 

Avi = Xv h Avi = Xvi. 

That is, the basis vectors of X no (X^) are the eigenvectors of A with zero in the ji-th, . . ., j p -th 
(zi-th, . . ., i m -th) components. □ 

In the rest of the paper we will denote the eigenvalues and eigenvectors for which (|2]) holds 
uncontrollable (respectively, unobservable) eigenvalues and eigenvectors. 

Remark 2.2 (Higher order integrators): The controllability (observability) analysis for k-th 
order integrators of the form x^ k \t) = aLx(t) + Bu(t) is equivalent to the analysis of the first 
order system ([!]). The statement follows, for example, by writing the conditions in Lemma |2T 



for the k-th order integrator system and exploiting its block structure. □ 

C. Motivating applications 

Next, we show three main areas of application for our results. 

Network of agents running average consensus: We consider a collection of agents labeled 
by a set of identifiers I = {1, . . . ,n}, where n G N is the number of agents. We assume that the 
agents communicate according to a time-invariant undirected communication graph G = (I, E), 
where E = {(i,j) & I x I \ i and j communicate}. The agents run a consensus algorithm based 
on a Laplacian control law (see e.g. H| for a survey). The dynamics of the agents evolve in 
continuous time (t G M>o) and are given by 

*<(*) = ~ 5Z W*) ~ s i(*))> » e {1, . . . , n}. 

For the controllability analysis, we consider a scenario in which some of the nodes have the 
possibility to apply an additional input that fully controls its dynamics. We call these nodes 
control nodes. This turns to be the model of a leader-follower network. For the observability 
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analysis, we imagine that an external processor (not running the consensus algorithm) collects 
information from some nodes in the network. We call these nodes observation nodes. In par- 
ticular, we assume that the external processor may read the state of each observation node. 
Equivalently, we can think of one or more observation nodes, running the consensus algorithm, 
that have to reconstruct the state of the network by processing only their own state. These two 
scenarios are captured by the model in equation ([T]). 

Remark 2.3 (Equivalence with other problem set-ups): Straightforward results from linear sys- 
tem theory can be used to prove that the controllability problem studied in flU and [9] and the 
dual observability problem studied in 031 can be equivalently formulated in our set up. □ 



Continuous time quantum and random walks: Dynamic systems induced by the Laplacian 
of a graph appear also in dealing with quantum and random walks [34]. We concentrate on the 
quantum counterpart of random walks, which have recently received great attention in the area 
of (quantum) information theory. The general idea of quantum information and computation is to 
solve common problems in information theory by using axioms and rules derived from quantum 
theory. Specifically, quantum walks are a computational variant of random walks in which the 
transition probability among the states follows quantum laws as opposed to standard stochastic 
laws. Formally, for quantum mechanical systems which are closed (i.e., not interacting with the 
environment) and finite dimensional, one considers the Schrodinger equation 

ihtp(t) = H(u(t))ip(t), 

where ip E C" is the quantum state and the Hamiltonian matrix H(u) is Hermitian and depends 
on a control u(t). Continuous time quantum walks are quantum systems whose dynamics is 
defined on a graph G. Specifically, the Hamiltonian has the form 



H(u) = H + Y^ ej-ej 



T 



J=H 



where H is the adjacency matrix or the Laplacian of a given graph G. The connection with our 
results appears for quantum walks on a grid (or lattice) graph with H being the grid Laplacian. 
The resulting dynamics is 



ihip{t) = (L + ^ / e j eju j (t))il>(t). 



j=n 
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In [|23l it is shown that the controllability of the above system, expressed by a Lie algebra rank 
condition, is equivalent to the controllability of the linear system Q. Our analysis is strictly 
related to the line pursued in ll23ll of finding more easily verifiable graph theoretic tests. 



Discretization of a class of partial differential equations: Next, we show how the discretiza- 
tion of partial differential equations containing the Laplace operator gives rise to an ordinary 
differential equation whose controllability and observability can be studied by using the tools 
developed in the paper. 

Let / : M. k — > M. be a twice differentiable real valued function, then the Laplace operator of 
/ is A/ := Yli=i da? " ^ n * s °P erator has a key importance in several physical phenomena. 
In particular, it appears in the heat and fluid flow diffusion, in wave propagation and quantum 
mechanics. Specifically, the density (temperature) fluctuations of diffusing material (heat) are 
described by the partial differential equation 

d ^^ = DA<f ) (x,t) + f(x,t), 

where <f>(x, t) is the density of the diffusing material (respectively the temperature) at location 
x G K fe and time t, D is the diffusion coefficient (respectively the thermal diffusivity) and f(x, t) 
is a material (heat) source. The wave propagation, arising in acoustics, electromagnetism and 
fluid dynamics, is described by the partial differential equation 



d 2 (j)(x,t) 



c'A(f>{x,t) + f{x,t), 



dt 2 
with 4>(x,t) the wave amplitude at x E IR fe and time t, c a constant, and f(x,t) a forcing term. 

If we consider a regular discretization of a d = 2 dimensional hyper-rectangular domain, e.g. 

[x^ in ,x^ ax ]x[x^ in ,x^ ax },for:i e {1, . . . , n x } and j e {1, . . . ,n 2 }, x [itj] = [x^ in +ih,x^ in +jh], 
where h is the discretization step, then the discretization of the Laplacian operator becomes 

where L is the Laplacian of an m x 712 grid graph and the vector (f) d (t) e M. ni ' n2 has components 
<j> d \ij]{t) = (/>(x[i t j],t). With this approximation in hands, the discretized versions of the partial 
differential equations above are ordinary differential equations with the same structure as in ([T|). 

III. Controllability and observability of simple eigenvalues in grid graphs 

In this section we characterize the controllability and observability properties of the simple 
eigenvalues of the grid, namely the eigenvalues of multiplicity one. 
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A. Laplacian eigenstructure of cartesian-product graphs 

An important property of graphs obtained as the cartesian product of other graphs is that the 
Laplacian can be obtained from the Laplacian of their constitutive graphs by using the Kronecker 
product of two matrices, see [28]. Given two matrices A G R. dxd and B G R. lxl , with [A]^ := ay, 
their Kronecker product A® B £ m dlxdl [ s defined as 

an-B dyiB . . . a\dB 

d2\B aiiB . . . a 2 dB 



A®B 

ddiB ayiB 
and their Kronecker sum as 

A@B = A 

Given the cartesian product of the graphs G\ 
Laplacian Lu of G\ □ . . . □ Gd is given by 

La — Li 



■ ■ a-ddB 

Ii + Id® B. 

. , Gd with Laplacian matrices L 1; . . . , L d , the 



L d . 



This structure on the Laplacian induces a structure also on its eigenvalues and eigenvectors. We 
state it in the next lemma, see ||28~1 . 

Lemma 3.1 (Laplacian eigenstructure of cartesian product graphs): Let G\, . . . ,G d be d G 
N undirected graphs and G = GiD . . . DG d their cartesian product. Let AJ, . . . , A£ be the 
Laplacian eigenvalues of the graphs G K and v*, . . . , v^ K the corresponding eigenvectors for k E 
{1, . . . , d}. The Laplacian eigenvalues and the corresponding eigenvectors of G are 



A, 1 , + A? + . . . + A" 



<2 



Id 



and 



< ® vl 



...®v, 



>d 



□ 



for «! G {1, . . . , m}, • • • , i d e {1, . . • , n d }. 

We are now ready to define a simple cartesian product graph. 

Definition 3.2 (Simple cartesian-product graphs): Let G and G' be two undirected graphs and 
let {Ai, . . . , Afc} and {A' l5 . . . , A' K } be the sets of distinct eigenvalues among all the Laplacian 
eigenvalues of respectively G and G'. We say that the graph Gn = G D G' is simple if the set 
{Ai + X' a | i G {1, . . . , k}, a G {1, . . . , k,}} contains only distinct eigenvalues. □ 

Using the associative property of the cartesian product the definition easily generalizes to the 
product of more than two graphs. 
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B. Controllability and observability of the simple eigenvalues 

We start with a lemma that relates the controllability (observability) of simple eigenvalues of 
a grid from a single node to the controllability (observability) of its constitutive paths. 

Lemma 3.3: Let Pi, ... , Pd, d E N, be path graphs of length respectively ni, . . . , rid and let 
G = PiD . . . nPd. Any simple eigenvalue A = Ai+. . .+\d of the grid graph G is not controllable 
(observable) from a node [(i)i, . . . , («)<*]» (*)i G {!,•••; ™i}> • • • » (*)d e {1? • • • > ^d}* if and only 
if there exists £ E {1, . . . , d} such that the eigenvalue \i of P? is not controllable (observable) 
from (i)e. 

Proof: From Lemma [27T the eigenvalue A of the grid graph is not controllable (observable) 



from a node [(i)i, ■ ■ ■ , (i)d] if and only if a Laplacian eigenvector w E R ni - n <* of A has zero 



[(i)i, . . . , (i)d] component. Using Lemma 3.1 and the assumption that A is a simple eigenvalue, 
any eigenvector w of A can be written as w — pv\ <g) . . . <g) Vd, where p G 1, p ^ 0, and Vi, . . . Vd 
are eigenvectors of the constitutive paths. Using the structure of the Kronecker product of d 
vectors, the [(i)i, ■ ■ ■ , (i)d\ component of w is zero if and only if at least one ve, £ 6 {1, ... , d}, 
has zero (i)i component. Indeed, {w)[(i) 1 ,...,{i) d \ = P( v i)h • • • • ■ { v d)i d which is zero if and only 



if there exists £ E {1, . . . , d} such that (ve) ie = 0. From Lemma 2.1 (ve) ie = if and only if the 



eigenvalue A^ of Pi is not controllable (observable) from (i)i, thus concluding the proof. 



From the previous lemma and Lemma 2.1 the next proposition follows straight. 

Proposition 3.4: A simple grid G = P\D . . . HP d is controllable (observable) from a node 
i = [(i)i, . . . , (i)d\ if and only if each path Pi is controllable (observable) from node (i)e- □ 

Using the property that any path graph is controllable (observable) from each external node, 
il5Tl . the next corollary follows. 

Corollary 3.5: Any simple eigenvalue of a grid graph is controllable (observable) from any 
corner node. □ 

We are now ready to characterize the controllability (observability) of simple grid graphs. 

Theorem 3.6 (Simple grid controllability (observability)): Let Pi, . . . , Pd be d path graphs of 
length respectively n\, . . . , rid and let G = P\\3 . . . DP^ be a simple grid. The following state- 
ments hold. 

(i) The grid graph G is not controllable (observable) from a node i = [(i)i, ■ ■ ■ , {i)d\, («)i £ 
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{1, . . . , rii}, . . . , (i) d G {1, . . . , n d }, if and only if there exists £ G {1, . . . , d} such that 



modp 



(^-(z),)^ P («)<-!) 



(3) 



for some odd prime p dividing n^; 
(ii) for any direction £ G {1, . . . , d} of G the following holds. For each odd prime factor p of ri£, 
the grid is not controllable (observable) from the set of nodes If = {i = [(i)i, . . . , (i)d] G 
I I (*)< = jp ~ 2 i i , J e {1, • • • , J}, and (i)i, . . . , (i)^_i, (i) /+1 , . . . , (i) d arbitrary}, with 
the following uncontrollable (unobservable) eigenvalues 



A,, 



7T 



2-2cos( (2i/ — 1)— ) +Ai + ... + A^i + A m + ... + A d , 



(4) 



and uncontrollable (unobservable) eigenvectors 



«V = «i <g) «2 ® •• ® ue-i ® Vui ® «£+!•• <8> Wd, K,£ as in (11) 



(5) 



with z/G {1, . . . , (p — l)/2}„ and A M , respectively w^, // 7^ £, any arbitrary eigenvalue, 
respectively eigenvector, of P M ; 
(iii) if a node i = [(i)i, • • • , (i)d] satisfies ([3]) for r < d distinct directions and, in each 
direction £, for kg < ne distinct prime factors, then the set of uncontrollable (unobservable) 
eigenvalues from node % is the union of eigenvalues with the following structure 



\i — A + Ai + . . . + Af_i + Xi + i + . . . + Arf, 
where each A is an uncontrollable (unobservable) eigenvalue of path Pi from (i) e and has 



the structure given in Theorem A.l (iv) and Remark A. 2 in Appendix. The corresponding 
uncontrollable (unobservable) eigenvectors can be built according to equation ([5]). □ 



Proof: The proof of statement (i) follows straight by combining the result of Lemma 3.3 



with the result in (i) of Theorem A.l 



To prove statement (ii), we start observing that, for each £ G {1, ... , off, the set of nodes 
if = {(i)ip — ^j^eii,..., 5 } is tne set of all nodes satisfying condition in Q for a given 



p in the path P t . Using Theorem A.l t we have that the uncontrollable (unobservable) eigen- 



values of Pi from this set of nodes have the form in ( [TO] ). Now, according to Lemma 3.1 
all the \ V £ as in Q are eigenvalues of the grid. Also, the corresponding eigenvectors are the 
ones given in (|5]). Using the result in Lemma 3.3 these eigenvectors have a zero in position 
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[0')i, • • • , (j)e-i, (i)i, (J)t+i ■ ■ ■ U)d] wit h {i)i satisfying © and (j) 1 , ..., (j)e-i, (j)e+i ■ ■ ■ (j) d 
arbitrary. To conclude the proof, we show that these are the only uncontrollable (unobservable) 
eigenvectors. To do that, we show that there exist nodes in Jf for which the only zero component 
is (i)i. For example, take any node with index [1, . . . , 1, (i)i, 1 ... 1] with {i)i satisfying Q. Since 
any path is controllable (observable) from the first node, the proof follows. 

Statement (iii) follows straight by combining the results in the previous two points with 
Theorem IA.1I ■ 

Next, we show, on the basis of the results in Theorem 3.6[ how to check the controllability 



(observability) of a simple grid from a given set of nodes or, equivalently, how to construct a 
set of control (observation) nodes such that the grid is controllable (observable). For the sake 
of clarity we present the procedure for a two dimensional grid (d = 2), however the procedure 
can be easily generalized to higher dimensions. 

First, we introduce some notation. Given two sets X = u£ =1 [xi :i ,, x 2>u ] and Y = U[, =1 [?/i iM , y 2 ^} 
with [xi iU , X2,v] 6 KxR and [yi >v ,y 2>l ,] G IR x IR, we say that XflF ^ if there exists 
[xi,d,x 2 ,u] e X and [3/1,^,1/2,^] E Y such that [x lfi ,x 2 , ] = [yi,p,, 2/2,/J, i-e. x ljP = y 1;P and 
%2,v = V2,p- 

Consider a two dimensional simple grid graph G = PiDi^ with P\ and P 2 of dimension 
ni and n 2 respectively. Let, for each I e {1)2}, n e = 2 m ' Yiif=iPl,v be a prime number 
factorization for some fc^GN and odd prime numbers p^i, . . . ,pi t k r Let I s = {ii, . . . , i m } be a 
set of control (observation) nodes with i a = [(i a )i, {ia)2\, a E {1, . . . ,m}. Now, we construct 
m sets Oi, . . . ,O m that will be used to define a simple rule for controllability (observability). 
For the sake of clarity we provide the rule to construct a set O a for a specific case. The general 
case can be easily deduced from the example. Suppose that i, satisfies condition ([3]) for pii and 
Pi j2 along direction 1 and for p 23 along direction 2. Now, define the set O a as follows. 

O a = [pi,i,P2,i] U . . . U [pi,i,p 2>k2 ] U [pi,2,P2,i] U . . . U [pi, 2 ,p 2 , k2 ] U [pi,i,P2;i] U . . . U [pi, kl ,P2,3\- 

We call 1; ... , O m a controllability (observability) partition of the set I s . 

The following proposition gives an easily implementable test for controllability (observability). 

Proposition 3.7 (Controllability (observability) test): Let G = P1QP2 be a simple grid graph 
and I s a set of control (observation) nodes with controllability (observability) partition 0\, . . . , O m . 
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Then G is controllable (observable) from I s if and only if 0\ n . . . fl O m = 0. □ 

Proof: The grid is controllable (observable) from the set I s if and only if the intersection 

of the sets of uncontrollable (unobservable) eigenvalues from each node is empty (equivalently 

if the intersection of the uncontrollable (unobservable) subspaces is the zero vector). Using 



statement (iii) of Theorem 3.6 the uncontrollable (unobservable) eigenvalues from each node 
are the ones in equation (|4j). Let A = Ai + A2, with Ai eigenvalue of P\ and A2 eigenvalue 
of P 2 , be a common eigenvalue to all the nodes in I s . Using the property in (iii), a node 
* = [W1W2] G Is can have that uncontrollable (unobservable) eigenvalue either because it is 
uncontrollable (unobservable) along one of the two paths or along both. Equivalently either Ai is 
an uncontrollable (unobservable) eigenvalue for P 1 from node (i) 1 and/or A 2 is an uncontrollable 
(unobservable) eigenvalue for P 2 from node (i) 2 . 



Now, using the result in statement (ii) of Theorem A.l all the nodes in a control (observation) 



set I \ are uncontrollable (unobservable) if and only if all of them belong to a set I\ = {£p — 
^Y~}e£{i,... 21 } f° r some factor p of n\ (ni being the length of Pi). The same argument holds for 
a control (observation) set I o2 on the path P 2 . This implies that the controllability (observability) 
on a path Pi can be studied by checking if all the nodes in the control (observation) set I oi share 
a common prime factor. 

Now, each point in the set O a , a G {1, . . . , m}, is of the form [pi,p2] where p\ and p 2 are 
prime factors of ni and n 2 (the lengths of P\ and P 2 ) respectively and at least one of the prime 
factors, say p x , is such that (i a )i G If 1 = {lp\ — ^S^lte-fi 211- Thus, each element in the set 
O a correspond to set of eigenvalues A Q = {X Ul ,u 2 G K I K lt u 2 = 2 — 2 cos ( (2^ — 1) — J + 
2-2cos((2z/ 2 -l)^),z/ 1 e {1,...,^} and v 2 G {1,...,^}} (with [% a ) x G If 1 or/and 
(i a ) 2 G If 2 ). The proof follows by observing that the controllability (observability) condition is 
that the intersection of the sets A Q be empty. ■ 

The following examples can be easily explained by using the proposition above. If at least one 
of the control (observation) nodes, say i\, is controllable (observable) in any direction, then the 
grid is controllable (observable). Indeed, the set 0\ will be empty. If all ni, i G {1, . . . , d}, are 
prime, then the grid is controllable (observable) if and only if one of the control (observation) 
nodes is controllable (observable) in any direction. Indeed, any O a , a G {1, ...,m}, can be 
either O a = {[ni, . . . , rid]} or O a = 0. 
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Next, we show a graphical interpretation of the controllability (observability) test based on 
the controllability (observability) partition. We present it through an example. In Figure [T] we 
show a two dimensional grid of length 7 x 15. It can be easily tested that this grid is simple. In 
each direction I G {1,2}, for each prime factor of n\ we associate a unique symbol to the rows 
(columns) of nodes that satisfy ([3]) for that prime number (in that direction). In particular, for 
the grid in Figure [TJ we associate a cross to the columns satisfying ([3]) for the prime factor 5 of 
15, a triangle to the columns satisfying ([3]) for the prime factor 3 of 15, and a pentagon to the 
unique row satisfying ([3]) for the prime number 7. 



o 




Fig. 1: Controllability (observability) partition for a 7 x 15 grid graph. 



Clearly, all the nodes that are not crossed by any line are controllable (observable). Also, 
a subset of nodes from which the graph is controllable (observable) can be easily constructed 
by suitably combining the different symbols. Equivalently, given a set of control (observation) 
nodes, testing the associated symbols easily gives the controllability (observability) property 
from the given set of nodes. For example, from the pair of nodes i\ = [1, 2] and %i = [4, 1] the 
grid is not controllable (observable). Indeed, the partitions are 0\ = [3, 7] and 2 = [3, 7] U [5, 7] 
whose intersection is [3,7]. The uncontrollable (unobservable) eigenvalues are, according to 
Theorem K^ 1 + (2 - 2 cos f ), 1 + (2 - 2 cos ^) and 1 + (2 - 2 cos ^). Following the same 
logic the grid is controllable (observable) form the set [1,2] and [1,3], but it is not controllable 
(observable) from [1, 2], [1, 8] and [4, 1] (and from any subset of them). We let the reader play 
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with it and have fun. 

IV. ElGENSTRUCTURE OF GENERAL GRID GRAPHS 

In order to characterize the controllability and observability of general grid graphs we need to 
exploit their eigenstructure. Indeed, the main difference with respect to the simple case analysis 
relies in the structure of the uncontrollable (unobservable) eigenvectors. While for simple grids 
they can be always written as the Kronecker product of two eigenvectors of the path (because 
the eigenvalues are all simple), this property does not hold for the eigenvectors of non- simple 
grids. Thus, the controllability (observability) analysis can not be performed by simply looking 
at how the zeros of the path eigenvectors propagate into the grid. Indeed, this analysis provides 
only necessary conditions for controllability (observability). 

This section will be organized as follows. First, we characterize symmetries in the structure of 
the grid eigenvectors. This analysis allows us to recognize the components of the eigenvectors that 
have to be equal. Second, we provide conditions to show what are all and only the components 
that are zero when a given component is forced to zero. Thus, with this results in hand, we are 
able to provide necessary and sufficient conditions for controllability (observability). 

We begin by characterizing symmetries of the path eigenvectors and then, using these results, 
we characterize symmetries of the grid eigenvectors by suitable grid partitions. For the sake of 
clarity we provide the analysis and results for two dimensional grids (d = 2). The results for 
higher dimensions are based on the same arguments and will be discussed in a remark. 

A. Symmetries of the path Laplacian eigenvectors 

We provide results on the structure and symmetries of the Laplacian eigenvectors of a path 
graph. The next lemma characterizes the symmetry of the path Laplacian eigenvectors. 

Lemma 4.1 (Symmetries of the path Laplacian eigenvectors): Any eigenvector v of the Lapla- 
cian of a path graph satisfies either v = Uv or v = —Uv, with II the usual permutation matrix. 
Proof: Let L e M nxn be the Laplacian of the path. Straightforward calculations show that 
L satisfies L = ilLil. Now, let v be a Laplacian eigenvector, then ilLilt; = Xv. Multiplying 
both sides by LT (and remembering that IT 2 = /), we get LUv = AIlt>, so that Uv is also an 
eigenvector of L associated to the eigenvalue A. Since any eigenvalue of L has multiplicity one, 
it must hold v = allv, for some nonzero a E R. Using the fact that the linear map II is an 
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isometry (i.e. it preserves the norm), ||iTu|| = ||i>||, it follows straight that either a = 1 or 
a = — 1, which concludes the proof. ■ 

In the rest of the paper we will denote S + (respectively S~) the set of vectors satisfying 
v = Uv (respectively v = —Hv). An important property of S + and S~ is that one is the 
orthogonal complement of the other, i.e. (S + ) ± = S~ . 

The next lemma relates the eigenstructure of a given path P to the eigenstructure of any path 
with length multiple of the length of P. 

Lemma 4.2 (Laplacian eigenstructure of P n and P kn ): Let Ai, . . . , X n be the eigenvalues of 
the Laplacian L n of a path P n of length n and vi,...,v n the corresponding eigenvectors. Then 
any path P kn of length kn, for some k G N, with Laplacian matrix L kn satisfies: 

(i) Ai, . . . , A n are eigenvalues of L kn ; 
(ii) each eigenvector Wi G M fcn of L kn associated to Aj, i e {1, . . . , n}, has the form 



111; 



Vf (UVi 



Proof: For the sake of clarity we prove the statement for k — 3, but the proof for the general 
case is easily generalizable. The Laplacian L kn can be written in terms of L n , whose structure 
is given in Appendix, as 

„t .. , / 



J kn 



J-'n * ^n&n 



_ T 



G\d n l^n T" Cl^i + CnC n Cn^l 







" e i e ;, 



T T 



T T T 

K K v i 



Lh/nWi 



Now, let us write the eigenvector Wi associated to \ as w^ 
v a in W l . Thus Wi satisfies 

L n v a + e n ((v a ) n - {v b )i) 

L n Vb + ei({vb)i - {v a )n) + e n ((v b ) n - (v c )i) 

L n v c + e 1 ((v c ) 1 - (v b ) n ) 
Now let us take v a = v i, v b = Hvi and v c — Vi, then 

e n ((vi) n - (rifj)i) 

ei((IlDj)i - (vi) n ) +e n ((Ilvi) n - (vi)i) 

ei((«i)i - (n« i ) n ) 



, with v a , v b and 



L kn Wi = XiWi + 



XiWi. 
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Last equality follows by the fact that (U.Vi)i = (vi) n and (U.Vi) n = (t>j)i (in general (Uvi 

(vi)n-e+i for any £ £ {1, ... , n}). 



Exploiting the result in the above lemma by using the result in Lemma 4.1[ it follows easily that 

T 



W: 



T T T 

vf vf vf 



for Vi = livi (and thus Wi = Hwi) and u>, ; 



T T T 

V f - V f v f 



T 



for Vi = — Uvi (and thus Wi = —Uwi). 



B. Symmetries of the grid eigenvectors 

Next, we provide tools to recognize symmetries in the grid eigenvectors, based on the graph 
structure, which will play a key role in the controllability (observability) analysis. 

Without loss of generality, let A = Ai,i + Ai,2 = . . . = A Mi i + A M) 2 be an eigenvalue of 
geometric multiplicity /i £ N, with Ai,i, . . . , A M! i (respectively Ai j2 , . . . , X^) eigenvalues of 
Pi (respectively P 2 ) and corresponding eigenvectors v 1 ,...v^ l (respectively W\, . . . , w^). The 
corresponding eigenspace V\ is given by 



m 



V x = {v £ R n ^\v = ^a^Vi ® Wi), on £ E}. (6) 

i=i 

As mentioned at the beginning of this section, it is worth noting that the eigenvectors in V\ 
do not necessarily have the structure of a Kronecker product of two eigenvectors (the set of 
vectors expressed as Kronecker product is not closed under linear combination). For this reason, 
in order identify all and only the zero components of these eigenvectors, we need to characterize 
their structure. 

Remark 4.3: For each node [is, £] such that the paths Pi and P 2 are controllable (observable) 
from v and £ respectively, all the basis eigenvectors of V\ have nonzero [u, £} component. □ 

Before stating the main results of this section, we need to introduce some useful notation. 
Given a path P n of length n £ N, for {i,j} C {l,...,n}, i < j, we denote P i: j the sub- 
path of P n with node set {i,...,j} (e.g., P 2:4 is the sub-path with node set {2,3,4}). Let 
G = Pi. ni nP m . n2 with P[. ni of dimension / • ni and P m .„ 2 of dimension m ■ n 2 . We call Gij = 
■P((»-i)ni+i):(*fn)n-P((j-i)na+i):yna). for »e{l,...,I} and j £ {1, . . . , m}, an m x n 2 sub-grid 
or a brick of G, see Figure |2j 

Let v £ M. l ' ni " m ' n2 be a vector of G, we call the sub-vector of v associated to G^ the vector 
Vij £ E nrn2 with components (%)[^], v £ {1, . . . , ni} and £ £ {1, . . . , n 2 }, given by (%)m = 

February 24, 2013 DRAFT 



19 



/// 



1 


Gn 


Gn 




(j~lm 












I 


Ga 


Gt2 




Glm 



Fig. 2: Partition of a grid into bricks 



(v )[(i-i) ni +v, (j-i)n 2 +e\- Informally, the sub- vector Vij off is constructed by selecting the components 
of v that fall into the brick G^ . 

Next, given a grid G = P ni DP n2 , with P ni and P„ 2 paths of length m and n 2 respectively, 
we introduce two useful operators that flip the components of a vector v associated to a grid G. 
Formally, consider the matrices 



^n ni ®/ n2 ; 



o 



712X712 



1 7l 2 



... 



7l 2 X n 2 



and (J ni ® II 



"-2 / 



n r 



o 



712X712 



Uri2X7l2 • • • -'■-'■712 



These operators flip respectively the first and the second sets of components. Formally, given 
a vector v G ]R nin2 associated to the grid G, with components (uW^, v E {l,...,ni} and 

£ G {1, . . . , n 2 }, let f ! = (Il ni eg) J n2 )f and t> 2 = (I m ® n„ 2 )t>. The vectors f i and f 2 are related 
to v by 

(vi)[ u ,e] = (v)[ ni - u+ ij], and (u 2 )[i/,<i = (")h« 2 -W]> 

for i/ G {1, . . . , rii} and £ E {1, . . . ,n 2 }. Finally, the composition of the two operators satisfies 
(II ni ®7 n2 )(7 ni ®Il n2 ) = (II ni ®n„ 2 ). Thus, when applied to a vector v, the composed operator 
flips both sets of components. That is, denoting v s = (Jl ni <g> H n2 )v, we have 

(V3)[u,(] = (l»)[m- I /+l,n 2 -^+l], 

for v e {!,... ,rii} and £ G {1, . . . ,n 2 }. 
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Lemma 4.4: Let G = P ni DP n2 with P m and P„ 2 paths of length respectively n x and n 2 . Any 
eigenvalue A of the Laplacian L of G is an eigenvalue of the Laplacian L of G = P\. ni DP m . n2 
for any I G N and m G N. 



Proof: From Lemma 3.1 each eigenvalue A of G (respectively A of G) can be written 



as A = A i + A 2 (A = A x + A 2 ) where A i (Ai) is an eigenvalue of P m (P/. ni ) and A 2 (A 2 ) of 



Pn 2 (Pm-n 2 )- From Lemma |4T2] all the eigenvalues of P ni (P„ 2 ) are eigenvalues of p. ni (P m . n2 ) 
so that the proof follows. ■ 

We are now ready to characterize the eigenvector symmetries by suitable brick partitions. 

Theorem 4.5 (Grid partition and eigenvector symmetries): Let Go — Pn 1 ^Pn 2 be a grid of 
dimension n% x n 2 with P ni and P n2 paths of dimension respectively rii and n 2 . Take any grid 
G = p. ni DP m . n2 of dimension Irii x mn 2 and let G^, i G {1, . . . , 1} and j G {1, . . . ,m}, be a 
partition into bricks of dimension ri\ x n 2 . 

Then for each eigenvalue (possibly non-simple) A of L : 

(i) A is an eigenvalue of L, and 

(ii) any eigenvector v of L associated to A can be decomposed into sub-vectors ify relative to 
the bricks Gij with 

Vij = (n ni ® i ri2 f~ l \i ni ® n^)^'- 1 )^ 

for i G {1, . . . , ni} and j G {1, . . . , n 2 }, where v is an eigenvector of L associated to A. 



Proof: Statement (i) follows straight by Lemma 4.4 To prove statement (ii), first, let us recall 
that the matrices (Tl ni ®I ri2 ) and (I ni £g>n n2 ) applied to the vectors v respectively flip the first and 
second components. Also, (II ni g) I n2 )^~ l ^w = w for i odd (respectively (I ril ®n n2 )( J_1 )w = w 
for j odd). Thus, we can just prove the result for (i, j) = (0, 1) and (i, j) = (1,0). 

Let /i be the geometric multiplicity of the eigenvalue A for the Laplacian L. Then, by 
Lemma |3.1[ a basis of the associated eigenspace V\ is given by \i vectors obtained as the 
Kronecker product of eigenvectors of the constitutive path graphs. That is, 



/' 



V x = {v G R ni - n2 \v = Y, a i( v i ® u>i), a i e R }- 



i=l 



where v { and W{, i G {l,...,/i} are eigenvectors of respectively P ni and P n2 associated to 
eigenvalues A^i and Aj2 such that A^i + \i )2 = A. Exploiting the Kronecker product and using 



the result in Lemma 4.2 we have 
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Clearly, the brick Gu coincides with the grid Go- Thus, we can compare the bricks Gij with 
the brick Gu. The sub-eigenvector corresponding to the "first row" of the brick Gu is given by 

On)[i,i] 



tUllJ[l,l:ni] 



adVihWi 



_(fn)[i, m ]_ 

Using the definition of brick components, the sub-eigenvectors corresponding to the "first row" 
of the grid Gi 2 is 



{Vl2)[l,lmi] 



(«ia)[i, 



l v 12j[l,m] 



fti^iITwi = fti(wj)i 



Wi 



(^i)i 



faOlLni] 



^llj[l,l] 



The proof for the other components follows exactly the same arguments. 

The above theorem has a nice and intuitive graphical interpretation, as shown in Figure [3] 
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Fig. 3: Graphical interpretation of Theorem 



4.5 
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Given a grid G and an eigenvector v associated to an eigenvalue A, we can associate a symbol 
to each node depending on the value of the eigenvector component. Then, we partition the grid 
G into bricks of dimension m X n 2 . Given the symbols in the brick Gn, the symbols in a brick 
G i: j, fori E {1, . . . , 1} and j E {2, . . . , m}, are obtained by a reflection of the brick Gij-i 
with respect to the horizontal axis, while the symbols in a brick Gij, for i E {2, . . . , 1} and 
j E {1, . . . , m}, are obtained by a reflection of the brick Gj-ij with respect of the vertical axis. 

Next, we analyze the eigenvector components of a brick whose dimensions are prime or, 
equivalently, the components of eigenvectors associated to eigenvalues that are not eigenvalues 
of smaller bricks. Recalling that any path eigenvector v satisfies either v = Uv (v E S + ) or 
v = —Uv (v E S~), we show how: (i) each basis eigenvector has a symmetry induced by 
the symmetry of the path eigenvectors, and (ii) the structure of a general grid eigenvector is 
influenced by the symmetry of the basis eigenvectors. 

Proposition 4.6: Let G = P ni DP n2 be a grid of dimension n\ x n 2 . For any eigenvalue A, let 
V\ be the associated eigenspace, with structure as in equation ([6]). Then, each basis eigenvector 
generating V\, (vi ® wi), satisfies one of the four relations: 

(v { ® Wi) = (n <g> i) (vt <s Wi) = (/ ® n) (vi <g> Wi) 
(vi ® Wj) = (n ® /)(^ <8> Wi) = -(/ ® n)(fi ® ^) 
(«< ® Wi) = -(n <8) /)(wi ® Wj) = (i ® n)(fi ® Wi) 

(^ ® Wi) = -(n (8) 7)^ (8) w,) = -(/ (g) n)(vj O Wi). 
Proof: Using the result in Lemma |4TJ v t and Wj belong either to S + or S~, that is, e.g., either 



Vi = Uvi or Vi = — Uvi. Suppose that, for example, Vi E S + and w, ; E S~ . Under this assumption 
[vi®Wi) = (rit>i(g>Wj) = (Uvi® Iwi) and, using the distributive property of the Kronecker product 
(vi®wi) = (U®I)(vi®Wi). Also, (t>i(g)Wi) = (fj(g)— ilwj) = {Ivi®— Uwj) = —{I ®W)(vi®Wi) . 
This gives the second of the four relations. The other three cases follow by the other three possible 
combinations of the symmetries of v* and Wi, thus concluding the proof. ■ 

In the following we denote the set of vectors satisfying each one of the four relations in the 
proposition respectively as S ++ , S + ~ , S~ + and S . A graphical representation of the four 
sets is given in Figure |4} We associate a symbol to each node depending on the value of the 
eigenvector component. Also, we denote with the same symbol but different colors nodes that 
have components of opposite sign. The result in Proposition 4^6] can be easily explained by using 
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(a) symmetry class S + ~ 



(b) symmetry class S~* 



(c) symmetry class S 




(d) symmetry class S 



Fig. 4: Graphical interpretation of Proposition 



4.6 



this graphical interpretation. Namely, each of the four cases in the proposition correspond to a 
scheme in Figure |4j 

The next remark gives an insight on the eigenvector components of the "central" nodes of a 
grid with odd dimensions. 

Remark 4.7 (Symmetries for grids with prime dimensions): If the grid has odd dimensions, 
ri\ and n 2 , the above symmetries have interesting implications for the nodes with components 
respectively ["^J], I e {l,...,ni}, and [u, Vf2 ^\, v e {1, ...,n 2 }. Indeed, the first set of 
components is zero for S~ + and S , while the second one is zero for S + ~ and S . □ 

This proposition has an important impact on the symmetries of general eigenvectors belonging 
to the same eigenspace (and in particular for each brick of a general grid). Clearly, if an 
eigenvalue is simple, then any associated eigenvector has the structure of a Kronecker product 
and thus one of the four symmetries. For a non-simple eigenvalue, any eigenvector of V\ can 
be written as a linear combination of the basis vectors, and thus, using the proposition, by the 
sum of at most four vectors each one having one of the four symmetries. Thus, in order to 
identify the symmetries of a general vector, we need to "suitably combine" nodes with the same 
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symbol and color in different classes. Three cases are possible: (i) if basis vectors of at least 
three different classes are present, by inspection in Figure [4] no symmetries are present, (ii) if 
all basis vectors belong to the same class, then also the linear combination does, and (iii) if the 
basis vectors belong to two of the four classes, a general eigenvector v (obtained as their linear 
combination) satisfies the following symmetries: 

a) (v)[ u ,£\ = ±(f)[„ 1 - J y + i i £] if the two classes share the first symbol (e.g., S ++ and S + ~); 

b) (v)[u,ei = ^ z ( v )[^n 2 -e+i] if the two classes share the second symbol (e.g., S ++ and S~ + ); 

c) (v)[v,e\ — =t(^)[ni-iH-i,n 2 -m] if tne tw0 classes share no symbol (e.g., S ++ and S ). 

A graphical representation of the above three symmetries is depicted in Figure [5] We associate 
the same symbol to nodes having the same absolute value of the eigenvector component. It is 
worth noting that we are interested in the absolute values because we want to classify all the 
components that can be zeroed simultaneously. 









— m)-(<> 







(a) first symbol in common (b) second symbol in common (c) no symbol in common 
Fig. 5: Symmetries of eigenvectors obtained as the sum of basis vectors of two different classes. 



V. Controllability and observability analysis of general grid graphs 

In this section we provide necessary and sufficient conditions to characterize all and only the 
nodes from which the network system is controllable (observable). First, we need a well known 
result in linear systems theory, see, e.g., [f33ll . We state it for the controllability property. 

Lemma 5.1: If a state matrix A E W nxn , n E N, has an eigenvalue with geometric multiplicity 
|i£N, then for any B E R nxm with m < fi the pair (A, B) is uncontrollable. □ 

The previous lemma applied to the grid Laplacian says that, in case the grid is non simple with 
maximum eigenvalue multiplicity fi, then the grid is not controllable (observable) from a set of 
control (observation) nodes of cardinality less than /i. 
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Using Lemma |2.1[ it follows straight that we can study the controllability (observability) 



properties of the grid separately for each eigenvalue. Namely, to guarantee controllability (ob- 
servability), we need to show that for each eigenvalue of the grid Laplacian L, there does not 
exist any eigenvector satisfying the condition in ([2]), i.e. having zero in some components. 

If A is simple, the corresponding eigenspace V\ in ([6]) is given by V\ — {v E IR 711 '™ 2 ^ = 
ai(fi <g> Wi),ai E IR}. Thus, finding the zeros of any eigenvector in V\ is equivalent to finding 
the zeros of the eigenvectors v 1 and u»i and propagate them according to the Kronecker product 
structure. Clearly, with this observation in hand, the analysis of any simple eigenvalue can be 



performed by using the tools provided in Section III 

For eigenvalues with multiplicity greater than one, next two considerations are important. First, 
not all the eigenvectors of A have the structure of a Kronecker product. Second, consistently 



with Lemma 5.1 it is always possible to find an eigenvector v G V\ with an arbitrary component 
equal to zero, for a suitable choice of the coefficients «j in (|6]). Thus, the controllability 
(observability) analysis does not depend only on the zero components of the path eigenvectors, 
but also on the symmetries in the grid eigenvector components. That is, for the eigenvalue under 
investigation, we want to answer to the following question. If we find an eigenvector with zero 
in an arbitrary component £, what are all and only the other components that are zero in the 
chosen eigenvector? We provide the analysis for non-simple eigenvalues of multiplicity two and 
discuss the generalization in a remark. 



On the basis of the eigenvector symmetries identified in Theorem 4.5 we can study the 
controllability (observability) of a brick. 

Next lemma provides useful properties of the eigenvector components in a brick. 

Lemma 5.2 (Polynomial structure of the eigenvector components): Let Go = P ni n\P n2 be a 
grid of dimension tii x n 2 . Then, any Laplacian eigenvector n = f®wof the grid, with v and w 
respectively eigenvectors of P n% and P n2 associated to eigenvalues A^ and X w , has components 
{u)[u,e\, v G {1, . . . ,ni} and I E {1, . . . ,n 2 } satisfying 

(i) {u)\y£} = p v (X v ) ■ Pe(X w ) ■ (v)i ■ (w)\, where p r (s) is the polynomial of degree (r — 1) 
defined as p 2 (s) = 1 — s for r = 2 and, denoting pi(s) = 1, by the recursion 

p r (s) = (2 - s)p r - 1 (s) - p r -2{s) (7) 

for r > 3; 
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(ii) if rii and n 2 prime, then p v (X v ) ^ and pe(X w ) ^ for any v E {1, . . . , rh ^} and 

Proof: First, notice that (u)\ Vt £\ = (v) u ■ (w)i. To prove statement (i), we need to prove that 
for a path of length m, any eigenvector v satisfies {v) v = p v {X v )(v)i for v E {1, ... , n{\. We 
prove the statement by induction. We exploit the eigenvector relation L ni v = X v v by using the 
structure of the path Laplacian given in Appendix. From the first row, it follows that (v )i — (v) 2 = 
\ v {v)\, so that (v) 2 = (1 — X v )(v)i. Then, from the rth row of the relation, we have —{y) r -\ + 
2(v) r — (f)r+i = X v (v) r and thus (v) r+ i = (2 — X v )(v) r — (i>) r -i- Plugging in the inductive 
assumption (v) r = p r {X v )(v)i and (v) r -i = p r -i(X v )(v)i, we have (v) r+ i = p r +i(X v )(v)i with 
Pr+i{X v ) = (2 — A„)pr(A„) — p r -i(X v ) which concludes the first part of the proof. 

Statement (ii) can be proven by showing that, for a path graph of length n\ with n\ prime, 
any eigenvector v has non zero components (v)i, . . ., (v) n t -i . This result is proven in fl7), thus 

2 

concluding the proof. ■ 

Next theorem gives necessary and sufficient conditions for two eigenvector components to be 
both zero in a brick. 

Theorem 5.3 (Simultaneous zeroing of eigenvector components): Let Go = P rai DP n2 be a grid 
of dimension rii x n 2 . Let A = A^i + Ai )2 = A 2 ,i + A 22 be an eigenvalue of multiplicity two, with 
Ai 5 i and A 2>i (A i)2 and A 2)2 ) eigenvalues of P ni (P m ). Let V\ be the associated eigenspace. Then 
there exists an eigenvector v E V\ with zero components [vi, £i] and [z/ 2 , £ 2 ], v\,v 2 E {1, . . . , m} 
and £i, £ 2 E {1, . . . , n 2 }, if and only if 

p Vl {Xi,i) -P4(Ai, 2 ) -Pv 2 (A 2 ,i) -pi 2 (X 2j2 ) = Pm(A 2 ,i) •Pi 1 {X 2 ,2) -^ 2 (Ai,i) • p^(Ai )2 ), (8) 

where p r {s) is the polynomial of degree r — 1 defined by the recursion in equation Q. 

Proof: To prove the statement, we look for an eigenvector v = ai(wi®wi)+ ot 2 (v 2 ®w 2 ) with 
ai and a 2 such that (f)^,^] = and (v)[„ 2t e 2 ] = 0- The condition (v)^/^ = is equivalent to 
<xi(vi) vi (wi)e 1 +a 2 {v 2 ) Ul {w 2 ) h = 0. From Lemma|i2|(i) we have ai^(Ai,i) • (vi)i -p^(A 2> i) • 



(itfi)i + a 2 p Ul (Xi^ 2 ) ■ {v 2 )i ■ pe 1 (X 2j2 ) ■ (w 2 )i = 0. Using the same calculations for the condition 
( V )W2M = 0' we can write the matrix equation 

Pu 1 (Xia)p£ 1 (X 1:2 )(v 1 ) 1 (wi) 1 Pu 1 (X 2tl )pe 1 (X 2:2 )(v 2 ) 1 (w 2 ) 1 
p V2 (Ai,i)p£ 2 (Ai, 2 ) (v i)i(wi)i p U2 (A 2 ,i)P4 (A 2)2 ) (v 2 )i(w 2 )i 



a i 

«2 



0. (9) 
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Since v is an eigenvector, a.\ and a 2 can not be zero simultaneously. Thus, the above equation 
is satisfied if and only if the matrix is singular. Imposing the condition that the determinant be 
zero gives equation ([8]). ■ 

Remark 5.4 (Extensions to higher grid dimension and eigenvalue multiplicity): If the grid is 
of dimension d > 2 the above theorem extends in a straightforward manner. Indeed, the result in 



Lemma pT2] (i) can be easily generalized as (u)[ iu ... lid ] = Pni^h) '■■■ 'Pitied) ' ( u i)i ' • • • ' ( v d)u 
with suitable adaptation of the notation, so that the condition in equation dHl) follows straight. 

If the grid has an eigenvalue of multiplicity /x > 2, then the theorem generalizes by con- 
sidering /i nodes. The condition in equation ([8]) follows by setting to zero the determinant of 
a /i x /i version of the matrix in ([9]) with elements (i,j) £ {1, . . . , //} x {1, . . . , //} given by 

P*<(Aj ) i)%(A il2 )(%)i(wi)i. D 

Next, we show a graphical interpretation of the controllability (observability) results obtained 



by combining the results of Theorem |3.6[ Theorem 4.5[ Proposition 4.6 and Theorem 5.3 We 



present it through an example. In Figure [6] we show a two dimensional grid of length 4x6. The 



analysis for the simple eigenvalues can be performed as explained in Section III This gives the 
cross symbol in the set of nodes [2, i] and [5, i], i 6 {1, . . . , 4} in Figure[6](b). Then, we partition 
the grid into bricks of dimensions 2x2 and 2x3. The eigenvalue Ai = 2 (respectively A 2 = 3) is 
an eigenvalue of multiplicity two in the brick 2 x 2 (2 x 3). The eigenvectors generating V\ 1 (V\ 2 ) 
belong to S + ~ and S~ + (S ++ and S ). This gives the symmetries in Figure [6] (a) according 



to Proposition |4.6| and the subsequent discussion. Using Theorem |5.3| it is easy to verify that all 
different symbols in Figure [6] (a) correspond, in fact, to distinct component values. Replicating 
the brick symbols according to Theorem 4.5 we get the structure in Figure [6] (b). Notice that 



in this particular case we have used the same cross symbol both for the non-simple eigenvalue 
A = 3 and for the simple eigenvalues. Finally, it can be easily tested that Ai = 2 and A 2 = 3 
are the only two non-simple eigenvalues. Given a set of control (observation) nodes, the grid is 
controllable (observable) if and only if the nodes do not have any symbol in common. If, for 
example, the control (observation) nodes share the top symbol, then the eigenvalue A = 2 (of 
the brick 2 x 2) is uncontrollable (unobservable). As for the simple case we let the reader play 
with the rule. 

To conclude, we provide a discussion on the importance and the effectiveness of the proposed 
methodologies in studying the controllability and observability (in general the dynamics) of 
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(a) bricks 2x2 & 2x3 




(b) brick partition 



Fig. 6: Graphical interpretation of the controllability and observability analysis. 



grid graph induced systems. First, we want to stress the fact that the proposed tools give 
strong insights on the structure and symmetries of the grid eigenvectors and of the controllable 
(unobservable) subspaces, as shown, e.g., in the above example. Furthermore, the proposed tools 
represent, clearly, an effective alternative to the standard tests in checking the controllability 



and observability properties. Notice that, the PBH test in Lemma 2.1 (to be performed for each 
eigenvalue) becomes prohibitive as the dimensions of the grid grow. Similarly, inspecting the rank 
of the controllability (observability) matrix is an operation that is ill-conditioned as the matrix 
dimension grows. As opposed to it, our tools involve the following operations. The test on the 



simple eigenvalues can be done simultaneously by using the tools in Section III and involves only 
arithmetic operations from number theory. The analysis for non-simple eigenvalues involves the 
following operations. First, using Theorem 4.5[ the grid can be partitioned into bricks of prime 
dimensions. This operation is based on a straightforward prime number factorization. Second, 
one has to compute the non simple eigenvalues for each brick. This can be done by using 



Lemma |3.1| and the closed form expression for the path eigenvalues given in Appendix. Third, 
for each multiple eigenvalue, one has to inspect the symmetries in each brick (again with simple 



operations on the brick dimensions according to Proposition 4.6) and the possible coincidence 



of symbols (by means of polynomial evaluations from Theorem 5.3). Finally, one should verify 
if there are multiple eigenvalues of the main grid that are not eigenvalues of smaller bricks. 
However, so far we have never found such a case in simulations, so that we conjecture that this 
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is unlike or even impossible to happen. 

VI. Conclusions 

In this paper we have characterized the controllability (by duality the observability) of linear 
time-invariant systems whose dynamics are induced by the Laplacian of a grid (or lattice) graph. 
We have shown that these systems arise in several fields of application as, in particular, distributed 
control and estimation, quantum computation and approximate solution of partial differential 
equations. We have characterized the eigenstructure of the grid Laplacian in terms of suitable 
graph decompositions and symmetries, and in terms of simple rules from number theory. Based 
on this analysis, we have shown what are all and only the uncontrollable (unobservable) set of 
nodes and provided simple routines to choose a set of control (observation) nodes that guarantee 
controllability (observability). 



Appendix 

In this section we briefly recall the results in [6], see also 0, on the controllability (observ- 
ability) of path graphs. The characterization of the controllability (observability) for grid graphs 
relies on these results. 

Since it is extensively used in the paper, we provide the expression of the path Laplacian, L n , 
and of its distinct eigenvalues Ai, . . . , X n 



1 


-1 







1 


2 


-1 









-1 


2 










-1 



and \ k = 2-2cos((k-l)-Yk = l,... 



n. 



The controllability (observability) of the path can be analyzed by using the PBH lemma 



in the form expressed in Lemma 2.1 First, it is known, H35I . that a path graph is always 
controllable (observable) from an external node (1 or n). Next theorem, which is Theorem 4.4 
in fl7), completely characterizes the controllability (observability) of a path by means of simple 
rules from number theory. 
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Theorem A.l (Path controllability and observability): Given a path graph of length n, let n = 
2 n ° Ylu=i Vv be a prime number factorization for some fceN and distinct (odd) prime numbers 
Pi, . . . ,pk- The following statements hold: 

(i) the path is not completely controllable (observable) from a node i E {2, . . . , n — 1} if and 
only if 

, .s modp , . , 

(n — l) = (t — 1) 

for some odd prime p dividing n; 
(ii) the path is not completely controllable (observable) from a set of nodes I s = {ii,...,i m } C 

{2, . . . , n — 1} if and only if 

,. ,i modp,. . >modp modp. . modp , . , 

for some odd prime p dividing n; 
(iii) for each odd prime factor p E {pi, . . . ,pk} of n, the path is not controllable (observable) 
from each set of nodes i| = {£p — ^y^gji... «} with the following uncontrollable 
(unobservable) eigenvalues 

\ v = 2 - 2 cos((2z/ -lAi/6{l r .., ^}; (10) 

and uncontrollable (unobservable) eigenvectors 



v u 



v T v o -(n^) T -^ o... (-i)HUv, 



(ii) 



where t>„ G IR^ 1 ^ 2 is the eigenvector of N( p _ 1 y 2 corresponding to the eigenvalue \ v for 
z/G{l,...,(p-l)/2};and 
(iv) if node i belongs to ij 3 = {£qj — ^5— }*£{i,...,— } f° r ^ — & distinct prime factors gi ^ . . . ^ 
5; of n, then the set of uncontrollable (unobservable) eigenvalues from node i is given by 

A. = 2-2co.(<2*-l)— ?— Y ^{l,..., ' 91 ---*)- 1 }. 
V 9i •••••«/ 2 

Also, the orthogonal complement to the controllable subspace, (AY)- 1 , (respectively the 

unobservable subspace, X no ) is spanned by all the corresponding eigenvectors of the form 



v v 



v T v -{Yiv v ) T -vl 0... (-i)?(n^) j 
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where v u E R^ 1 '-'*) i)/2 j s ^ eigenvector of N^ qv ,__. qi yxy 2 corresponding to the eigen- 
value A„ for v E {1, ... , ((<?i -...-«)- 1 )/ 2 )- 

D 



Remark A.2 (General version of Theorem A.l): In the general case of a path graph of length 
n = 2 n ° Yl u=1 Pu, where px,...,pk are not all distinct, statement (i) and (ii) of Theorem A.l 
continue to hold in the same form. As regards statement (iii), it still holds in the same form, 
but it can also be strengthen with a slight modification. That is, for each multiple factor p with 
multiplicity k, the statement continues to hold if p is replaced by p a with a E {l,...,k}. 
Statement (iv) holds if for each prime factor p with multiplicity k we check if node i belongs 
not only to If but also to each If with a E {l,...,k}. Consistently the uncontrollable 
(unobservable) eigenvalues and eigenvectors considered in the statement must be constructed 
by using p a instead of p, where a = max a {a E {1, . . . , k}\i E if}. □ 
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